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1. Introduction 

The BFKL Pomeron calculus [|, @ gives the simplest and the most transparent approach to the high energy 
amplitude in QCD. This approach is based on the BFKL Pomeron || and the interactions between the 
BFKL Pomerons (§ |, |, |], which are taken into account in the spirit of the Gribov Reggeon Calculus 
||. It can be written in the form of the functional integral (see Ref. |J) and formulated in terms of the 



generating functional (see Refs. M, [LOL O, 12, O]) 



In Ref. [14] 1 (see also Ref. [|T5|) it is shown that the BFKL Pomeron calculus can be reduced to the 



following functional integral 

Z[$,$+] = J D<Z>{b,k;Y)D<t> + {b,k;Y)e s[ ' s ''* +] (1.1) 



with 



S = So + S! = J d 2 bd 2 k$+(b,k;Y) ( ^hpll _ g. (J d 2 k'K(k,k')<S>{b,k';Y) + (1.2) 



^ee Eq.(3.18), Eq.(3.55) , Eq.(4.97) and Eq.(4.98) of Ref. @ for more details. 



- 1 - 



$+(&, k; Y) $(6, k;Y) + a 2 s {$ (6, fc; F) (6, fc; F) - $ + (6, fc; F) $ (6, fc; F) $ (b, fc;F)})) 
where X (&, fc') is the BFKL kernel in the momentum representation, namely, 

<D(6,fc;F) 



y d 2 k' K(k,k') $(b,k';Y) = j d 2 k 



(k - k'f 



(k - k') 2 {k n + (Jfe - /c') 2 ) 



1.3) 



The theory determined by the functional integral of Eq. can be written [14] in the equivalent 

form introducing the generating functional Z(Y, [u(b,k)]) defined as M 



Z(Y - F ; [u(bi,ki)}) 



;i-4) 



/n 
P n (F - F ; 6i, fcj; . . . ; 6j, fcj; . . . ; 6 n , fc„) u(6», fcj) d 2 ii 2 fc, 

ra=i i=l 

where n(6j, /cj) are arbitrary functions and P n (F — Yq; b\, k\; . . . ; 6j, kf, . . . ; 6 n , /c n ) is the probability to find 
n-dipoles with transverse momenta ki- 



For this functional we have the following evolution equation [14] 
dZ(Y-Y ;[u(b,k)}) 



dY 



X [u(b,k)] Z(Y -Y ;[u(b,k)}) 



with 



xm 



d 2 bd 2 k f ^ ( - / d 2 k' K{k, k')u(b, k') - £ , - + u(b, k) u(b, k) " 



6u(b, k>) 



8u(b, k) 



1 



-q s (u(b,k) u(b,k) — u(b,k)) — - 



2 5u{b,k)5u{b,k) J ' 



(1.5) 

(1.6) 
(1.7) 



In this paper our goal is to solve Eq. ( |1.5| ) deeply in the saturation region or, in other words, for the 
dipoles which sizes are much larger than the typical saturation size 1/Q s (x). In the next section we set 
the problem and will discuss the main assumptions that simplify the problem. In section 3 we discuss the 
mean field approximation. In this approximation the expression for x[ u ] 1S simple and it is determined by 
Eq. (|1.6| ). Solution in this approximation is known (see Ref. [ji7|] ) but we will solve this equation using a 
different method. Therefore, we view this section as a training ground for our method of finding a solution. 

In section 4 we focus our efforts on searching a solution to the general equation (see Eq. (|1.5|)). We 
consider the solution which shows the geometrical scaling behaviour in the saturation domain and which 
matches with the linear evolution equation at r ~ 1/Q s (x) where r is the dipole size and Q s (x) is the new 
scale: saturation momentum. 

In conclusions we summarize the main results of the paper and compare them with other attempts to 
solve this equation . 
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2. Problem setting and the main assumptions. 

The problem that we would like to solve, is the scattering of very small colourless dipole , which size is 
r, with the large target (let say with a nucleus with radius R, R 3> r). We are interested in finding the 
asymptotic behaviour of this amplitude deeply in the saturation region where «5ln(r 2 Q 2 ) 1. 

Generally speaking , we need to fix initial and boundary conditions to solve Eq. (|1.5| ). These conditions 
are well known and they have the form 



Initial conditions: at Y = Z (Y = 0; [u(b, k)]) = u(b, k) ; 
Boundary conditions: at Ui = 1 Z (Y; [u(b, k) = 1]) = 1 ; 



(2.1) 
(2.2) 



Eq. ( |2.1D means that at low energy the interaction can be reduced to the single BFKL Pomeron exchange 
while Eq. ( |2.2[ ) follows directly from the meaning of P n as probability to find ra-dipoles and reflects the 
normalization for the probability (5^-P ra = !)• 

In terms of anomalous dimension 7 we replace the general BFKL kernel (see solid line in Fig. [I] 



if(M')e (7 - 1)hl(fc ' 2/fe2) = X(7) = 2^(1) - V(7) - ^(1-7) 



(2.3) 



by 



xin) 



1-7 



for r 2 Q 2 S < 1 ; 
for r 2 Q 2 > 1 ; 



(2.4) 



(see dotted-dashed line in Fig. [T]) 

This approach is quite different from the diffusion approximation for \ (7) (see dotted in Fig. [I] 

X in) = 21n2 + 28C(3)( 7 -^) 2 



(2.5) 



^3], |2^] to reduce the general equation to so called Fisher-Kolmogorov- 
(F-KPP equation). (,{x) in Eq. (|2.5[) is the Riman zeta function. 



which was used in Refs. 15, 22, 
Pertrovsky-Piscounov equation 

The kernel of Eq. (|2.4j ) sums the contributions of the order of (as ln(r 2 A 2 ))™ for r 2 Q 2 <C 1 where 
as ln(r 2 A 2 ) 3> 1; and it leads to summation of the terms of the order of (as ln(r 2 Q 2 )) n in the kinematic 
region where r 2 Q 2 S » 1 and as ln(r 2 Q 2 ) » 1 |l^, f7|. 

It means that this kernel can describe the behaviour of our systen deeply in the saturation region 
(r 2 Q 2 S 3> 1) or/and in the perturbative QCD region at r 2 Q 2 -C 1 while we cannot expect that our 
approximation will be able to reproduce the exact solution in the vicinity of the saturation scale r 2 Q 2 ~ 
1). These features of this kernel one can see from Fig. [l|. Fig. [l] shows that the approximate kernel of 
Eq. Q2.4p can describe the exact BFKL kernel only at 7 close to and 1 while it leads to a contribution 
which considerably differs from the exact kernel at 7 = 1/2. For 7 close to 1/2 we can use the diffusion 
approximation of Eq. (|2.5| ) which fails to reproduce the correct behaviour inside the saturation region. One 
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can see this since we cannot obtain the correct asymptotic behaviour |17J in the mean field approximation 
with the kernel of Eq. Q2.5p . On the other hand, the kernel of Eq. ( |2.4| ) cannot describe the behaviour 
of the scattering amplitude in the region close to the saturation scale. In particular, Eq. ( [2.4| ) gives the 
saturation scale as 



HQ 2 S R 2 



4a s Y = 4a s ln(l/x) 



where as = (N c /ir) as, while the correct saturation scale is equal [Tj| [0| 



ln{Q 2 s R 2 



X (icr) 
1 - 7cr 



with 7 cr determined by equation: 



X(lc 



1 - 7c 



djcr 



Eq. gives 7 C 

(see Fig. |). 



1/2 that reduces Eq. (2.7) to Eq. (2.6) while the correct value of j C r is 7, 



(2.6) 

(2.7) 
0.37 




We used two additional assumptions to reduce the 
general functional integral for the BFKL Pomeron cal- 
culus given in Ref. H to the simple form of Eq. (|l.2| ). 
Namely, we assume that at high energy we can use the 
BFKL Pomeron at 7 close to 1/2 which leads to 



L 



12 



T 4 
x 12 
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Ycr Y 1 "Y cr 

Figure 1: The kernel of the BFKL equation: the 



solid line shows the exact expression of Eq. (2.3); the 
kernel of Eq. (2.4) is plotted in dotted-dashed line; 
and the diffusion approximation of Eq. ( |2.5| ) for the 
BFKL kernel, that leads to F-KPP equation, is pre- 



Vy G(x 1 ,x 2 -Y\x' l ,x' 2 ;Y') = 

= G{x 1 ,x 2 ;Y\x' 1 ,x' 2 ,Y') (2.8) 

where G (x\,x 2 ;Y\ ; Y') is the Pomeron Green's 
function. 

This assumption works quite well both for diffusion 
approximation for the BFKL kernel (see Eq. fl2.5| )) and 
for the kernel of Eq. ( |2.4| ) since 7 cr for this kernel is 
equal to 1/2 (-y cr = 1/2). 



It should be mentioned that Eq. (1.2) is written 
assuming that 



sented by dotted line. 7 cr is calculated using Eq. (2.7) 



X\ + x 2 



> xi - x 2 ; 



(2.9) 



Indeed, this assumption looks natural for the mean field approximation, since we expect that the typical 
size of the dipoles will be of the order of 1/Q s (x) (Q s (x) is the saturation scale) while the typical impact 
parameter of the scattering dipole should be much larger (at least of the order of the size of a target (i?) 
). However, in the Pomeron loop we have integration over b and we should study the b dependence more 
carefully. Generally speaking, the integration over impact parameter in the Pomeron loop leads to the value 
of b of the order the typical size of the dipoles in the loop and can change considerable the conclusions of 
the approaches where the b dependence has not been taken into account (see p0, 22, S3, 153]). 
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For our kernel, that sums log contributions, the impact parameter dependence does not influence 
on such summation since the BFKL Pomeron shows the rapid fall down as a function b (oc 1/6 4 ) and 
integration over b does not generate any logarithmic contribution. The b dependence for our kernel has 
been discussed in details in Refs. [I], |17| and we will show how it works in the next section. 

It should be mention that the relation between fields in momentum and coordinate representation 
looks as follows 

<f>(xi,x 2 ;Y) = x\ 2 J d 2 ke itSl2 <S>(k,b;Y) (2.10) 

where x±2 = x\ — xi- 

We would like to solve directly Eq. (|1.5| ). Our main assumption is that Z (6, k; Y) actually is a function 
of one variable Z (£; Y, [u] ) 2 : 

e = ln(Q2(z)/fc 2 ) (2.11) 



In finding the solution we replace the generating functional of Eq. (1.4) by the generating function 



Z(£,u) = Y, p n(0u n (2.12) 

n=l 

We can use Eq. ( |2.12| ) instead of the generating functional of Eq. ( |1.4[ ) since the scattering amplitude 



is determined by the following equation pO, 11 1 



N(Y;k;[ 7i }) = - £ i-i- ln ({h}; b; Y ; R) ]J d 2 h —Z (Y - Y ; {ki}, u) \ u - 



n=l 
oo 



?»1 

- - 7 



V ^-^ ln ({k i = l/R};b;Y ;R) Z(Y-Yo;{ki = l/R 2 }, U )\ u=1 = 

n=l 

00 ( — \\ n 

V j n ({ki = 1/R}; b; Y ; R) Z (f; u) \ u=1 (2.13) 



n 

n=l 



7n({ki}', b;Yo;R) in Eq. ( 2.13| ) is the scattering amplitude of n dipoles with momenta ki with the target of the 



size R at low energy. In logarithmic approximation we can consider this amplitude as 7 n ({fci}; b; Yq; R) = 
7n({fc« = b; Yq; R) fl™ S^ 2 '(ki — 1/R)- Therefore we need to know only generating function of 

Eq. ( |2.12| ) to find the scattering amplitude. 

In the perturbative QCD region (r 2 Q 2 (x) ~ Q 2 /k' 2 <C 1) the scattering amplitude has the following 
form for the kernel of Eq. ( |2.4D 



N(Y = ln(l/x);ln(k' 2 R 2 )) = a| exp ( 2 y/a s Y ln(k' 2 R 2 ) - \n(k' 2 R 2 ) ) (2.14) 
which can be translated in the initial condition, namely, 



2 For simplicity we neglect the dependence of Z on impact parameters. For large nuclei this is a good approximation. 
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Z (Y;ln(k' 2 R 2 ); u) =1 - (1 - u) exp ( 2 y/a s Y ln(k' 2 B 2 ) - \n(k' 2 R 2 ) 



(2.15) 



In Ref. [26] it is proved that the solution of the linear evolution equation in the vicinity of k 2 — > Q 2 shows 
the geometrical scaling behaviour and in terms of the generating function this contribution looks as 



Z(£; U ) 



k<o, i?i<i 



(1 -u) exp ( - £ 



(2.16) 



for the solution of Eq. ( |2.15| ). In Ref. [26] it is shown that the dependence of the perturbative solution on 
the only one variable £ is a general property also for the exact kernel of the BFKL equation. Eq. fl2.16| ) 
will be used as the initial condition in our solution. 



To check the strategy of our approach, especially Eq. ( |2.12| ) , we first consider the BFKL Pomeron 
interaction in the kinematic region near to the saturation scale where we expect the geometrical scaling 
solution of Eq. ( |2.16| ). We will show that in leading log approximation of perturbative QCD which corre- 
sponds to the kernel of Eq. (p~4|), we can avoid the assumptions given by Eq. (2.8) and Eq. (2.9), while 
Eq. ( 2.12| ) will be proved. 



3. The BFKL Calculus in the geometrical scaling region 

We start with the simplest enhanced diagram of Fig. ^-a. We will write the expression for this diagram for 
function <t> which is defined as 

r d 2 b e 1 ^ N(^,p;b) (3.1) 



<t>(£;p,q) = 

where £ is defined by Eq. (2.11) and p = ln(k 2 R 2 ). The diagram can be written as follows 



<MFig.@-a|£,p;<z) = J d?d p'd 2 q' d£" d p" <j> (£-£', p - p' , q) 

G 3 p , 4> (£' - £" , pi - P" , g- if) <P (£' - £" , P ' - P" , <f) G 3P 4> (£" - £ , p" - Po, q) 
In Eq. (|3.2[) we introduced the variables £ instead of rapidity y and p instead of k 2 . 



(3.2) 



The integration over q' in Eq. (3.2) has been discussed in details in Refs. ||, £7]]. The result of this 
discussion is the fact that the integration over q' sits at q' = the lowest virtuality in the loop. In other 
words, f d 2 q' = irk" 2 in Fig. ^-a if we calculate this diagram in the region of £ < 0. Taking this into 
account we reduce Eq. ( |3.2| ) to the following one 

0(Fig.@-aK,tf<z = O) = ^ f J" dtfdp' J* j' d?>dp» <£(£-£', p-p J ,q = Q) (3.3) 

= 0) </>(£' -£'V-/A<? = 0) 0(£"-£o,p"-/oo,<? = O) 
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It is clear the summation all diagrams can be done using the following functional integral 

Z[$,$+] = J D^,p)D^ + (C,p)e S[ ^ +] 



(3.4) 



with 



S = S + Si = f d£ dp $+(£,/,) U 8 ®f^ P) ~ f 



dp' 



4^1 + Htp ') + (3.5) 



+ 4>+ (5, p')* (tJ) +a' s {9{(,fl) 9 + {(,//) - 3>+ {(,/) * (£,/)})) 

As was shown in Ref.[|l4|] this functional integral can be rewritten as the equation for the generating 
functional , defined as follows 



namely, 



dp / dp 



/n 
P n (£, pi; ■ ■ ■ ; pr, ■ ■ ■ ; pn) J\ u (p^ d P i 

71=1 1=1 

az(£,M) 



(3.6) 



(3.7) 



(u(p) - u{p) u(p') 



a s (u(p)u(p) -u(p)) 



s 2 z(t [«]) 



This functional equation gives the full description in double log approximation for £ < 0. In particular, 
it leads to the linear equation for the functional N = 1 — Z (£, u = 1 — 7) if we neglect the terms of the 
order N 2 . This equation is 



d 2 N dN 
+ 



d£dp di 

It is easy to check that N of Eq. ( 2.14j ) satisfies this equation. 



N 



(3. 



However, Eq. (3.7) can be simplified in the region of |£| 3> 1 and p 3> 1 in the kinematic region where 
we have the geometrical scaling solution of Eq. ( p. 16 ). In this region Eq. fl3.3| ) has a simpler form, namely, 



■ (Fig. @-a|£) 



n 7 



dp' I dc I d^ 0(e-o ^(e'-r) (3.* 







^5 

16 





[ cdedc [ d^ #(£-o 0(£'-r) ^(e'-r) <H£"-&) (3.10) 

Co JO 



In calculating Eq. ([?!]) we restrict ourselves by the kinematic region £ < where we can safely use 
Eq. (12.160 for the BFKL Pomeron. All integrations in Eq. (Q and Eq. fl3.10|) are taken in the so called 
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double log approximation which corresponds to Eq. ( |2.4| ) for the kernel and Eq. ( |2.15D for the solution of 
the linear equation. We also used the kinematic restriction p < £ that follows from Fig. Q-b. 

Eq. ( p. 10 ) suggests that the problem of summing all BFKL Pomeron diagrams in the saturation region 
for £ 3> 1 can be reduced to the solution of the following equation for the generating function defined by 
Eq. ( gig ) and Eq. flO| 



dZ&u 



■ £ u (1 — u) 



dZ 



+ a s £ u (1 — u) 



d£ " du 

This equation can be rewritten in more convenient form introducing y = £ 2 / 



du 2 

I, namely, 



dZ (yyu] 

dy 



u(l 



u 



dz{y-u\ 

du 



+ a 2 s u{l 



u 



d 2 z{y-u) 

du 2 



(3.11) 



(3.12) 



Eq. ( 3.12[ ) is the same as for the BFKL Pomeron calculus in zero transverse dimension (see Ref. [29] 
and references therein). This equation has been solved [g^] . We will discuss below the main propertiers of 
the solution to Eq. (|3.12j ) but first we start from the solution to our problem in mean field approximation 
to check the main assumption that we make on the way of obtaining Eq. ( [3.12 ). 



4. Exact solution and the mean field approximation (MFA) 

In the mean field approximation x[ u \ is determined only by Eq. ( |l.6|) . This fact considerably simplifies the 
equation since such an equation belongs to the well known class of the Louisville equations. The Louisville 
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equation has a general solution, namely, Z (Y — F(u)). Using this form for the general solution as well as 
two conditions (initial and boundary, see Eq. ( |2.1| ) and Eq. ( |2.2| )): we can easily prove that for Z we have 
the non- linear equation M, ffO[, namely 



dZ(Y-Y ;b,k] 
dY 



g (- J K(k, k') Z(Y- y ; b,k') + Z 2 (Y - Y ; b, k )) (4.1) 



This equation leads to Balitsky-Kovchegov equation [£(], 21] for the scattering amplitude and it has been 
solved in Ref. [17| and the solution turns out to be the function of the only one variable ^ (see Eq. ( 2.11| )). 
In Ref. [17] it is shown that at large £ the generating functional Z (Y — Yq; 6, k! ) for the kernel of Eq. ( |2.4[ ) 
behaves as 



Z(Y-Y ;b,k') = Cexp(-^pj 



(4.2) 



In the log approximation for the BFKL kernel in the saturation region which is related to the kernel of 



Eq. ( |2.4D , (see Ref. [17] and references therein) we have : 



r r Q2 ^ dk 2 

J d 2 k' K (k, k') Z (b, k'; Y) J ^-Z(b,k';Y) 



(4.3) 



In the MFA Eq. (3.12) reduces to very simple Louisville equation, namely 



dZ (y-u 



dy 



u(l — u) 



dZ (y;u] 
du 



which has the solution 



z(y-,u, 



1 + (1 - u) (exp (y) - 1) 



(4.4) 



(4.5) 



One can see that Eq. (4^) satisfies the initial condition of Eq. ( 2.16 ) if we define 3^ = £ 2 /8 — £/2. 
For £ — > the scattering amplitude is equal to 



N(y, 1 = i-u) 



7 exp (y) 



1 + 7 (exp (y) - 1) 



(4.6) 



Since 7 is of the order of as we obtain Eq. ( p. 14 ) at small £ with our definition of y . Of course, we cannot 
guarantee term ^£ since we assumed that £ 3> 1, therefore, we add this term to guarantee the matching 
with the initial conditions at ^ — > 0. 



Therefore, we showed that our procedure involving Eq. (2.12) and Eq. ( 2.16j ) leads to the correct 
behaviour of the amplitude in the mean field approximation. Therefore, we confirm our main assumption 
that for searching for behaviour at large values of £, we can consider of generating function of Eq. ( p. 12 ) 
instead of the generating functional of Eq. ( |1.4[ ) . 
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5. A possible solution with the geometrical scaling behaviour 



5.1 General properties 

Eq. ( 3.12Q has been solved in Ref. [29|., where it shown that the solution for the scattering amplitude 
approaches at high energy behaves as 

AT(£) ^ i _ exp(-C(K)e 2 /2) (5.1) 



Ref. where k = l/ar|. The exact value of the factor C(k) is found in Ref. [29], however, here we will 
illustrate the main property of the exact solution using an approximate semi-classical method developed 
in Ref. which explores the large value of k. 

First we find the asymptotic solution to the master equation from the following equation: 

, dZ(^u) 1 , , d 2 Z (£;u) . N 

The solution that satisfies the boundary condition of Eq. ( fE^ ) is equal to 

I — e KU 

Z asymp (y = oo,u) = 1 _ eK (5-3) 

The second step is to search the solution in the form 

Z (£; u) = Z asymp (y = oo, u) AZ (£; u) (5.4) 

assuming that d 2 AZ (y; u) /du 2 <C k (, dZ asymp (y = oo, u) /du ) (dAZ (y-,u) /du). Since k » 1 this 
assumption looks reasonable. For AZ (£; u) we obtain the equation 

^ (y u ) dAzjyu) dz asymp (y = ^,u) dAz (y u ) 

dy ~ ~ u ^- u ) ^ + 2u ^- u > ^ du (5 ' 5) 

Assuming that k 3> 1 Eq. (|5.5[) reduces to the form 

dAz (y-u) ^ dAz iy-u) 

dy = u{l ~ u) — du — (5 ' 6) 

Eq. ( |5.6| ) has the solution 

AZ(y-u) = $ (y + Hjz^) (5 - 7) 

with function that should be found from the initial conditions of Eq. (|2.1| ). 
The solution is 

1 — e KU u 

Zasymp (V = oo; u) AZ (3?; u) = u- -~ with u = — — — (5.8) 

y F 1 — e KU u + (1 — u) exp (—y) 
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However we cannot use directly Eq. (|5.8| ) to determine the asymptotic behaviour of the scattering amplitude 
at high 3^ since we need to find the typical values of u for such an amplitude at large values of £ or to 
find the solution in the entire kinematic region. The latter we cannot do since our equation is valid only 
at £ 3> 1. To solve this problem at £ 3> 1 we need to use the i=channel unitarity constraints ( see Refs. 



p3| , 34, |13|). we need to use the t-channel unitarity constraint ( see Refs. |33], 34, |13[), which has the 
following form 

(0 - t ^ 7" (5.9) 

n=l 

where j BA is the amplitude in the Born approximation for the interaction of two dipoles at low energy, 
which is equal to a| = 1/k in our case. Using k 3> 1 we see that in this limit ^§ = (K^) n Z and 
£ = exp (-y) at u = 1. 

Using these observation we obtain that 



i^zVgi, i _ e -4« 2 _ 2k e -? 2 /32 (5 . 10) 



The choice £' = £/2 is clear since it gives the slowest approach to N = 1. Eq. ( 5.10 ) can only be considered 



as a approximation to the exact solution of Ref.[29] . It should be stressed that Eq. ( |5.10 ) leads to quite 



a diffrent result in comparison with the Mueller-Patel-Salam-Iancu approach [33]. However, the MPSI 
approach has an advantage to be correct also for £ ~ 1. This is the reason why we consider our problem 
in this approach in the next section. 

5.2 The Mueller-Patel-Salam-Iancu (MPSI) approach 



This approximation has a legitimate region of applicability 34]; can be used for small values of £ and 
it describes quite well in the entire kinematic region the exact solution to the BFKL Pomeron calculus in 
zero transverse dimension 



In this approach the scattering amplitude can be calculated using the following formula |34| 

N MPSI (x,y;x',y'-Y) = (5.11) 

j ^ 6 n Z MFA iY _ Y ,. [ui])i 5 n Z M FA{Y ,. [vi]) BA 

\Vi=l 1 \X%, t/j, Xj, IJi) 



^ (-l) n n! Y\ dxidyidx'idyi 
n=l ^ i=l 



where j BA (xi, yf, x[, y[) is the amplitude of interaction of two dipoles in the Born approximation of per- 
turbative QCD. 

In the case of the geometrical scaling solution the MPSI approximation formula has an elegant form 
[j36| , namely 

N MPSI (6 = 1 _ | exp (-^ A. £\ N MFA (e _ e . u) N MFA ^ y) 1 (fU2) 
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where the generating function j\r (£> u ) is defined as 



N 



MFA 



1 



(5.13) 



The generating functional in MFA for the kernel of Eq. ( |2.4f) has been found in Ref. [17|. In this paper 
it is shown that the solution to Eq. ( f4.l|) has the form that in our notations looks as follows 

AT(£,7=l-«) = 1 - exp(-C(£,u)) 

where function C(£i u ) can be calculated from the implicit equation 



e = V2 

For small values of £ — > 7 Eq. (5.15) gives 



+ (exp(-C') - 1) 



(5.14) 
(5.15) 



ln(C(e,«)/7) 



-£ which leads to ((£,,u) 



7e2 > 



7(1 + 5 



(5.16) 



It is easy to see that Eq. (|5.14Q with £(£,«), determined by Eq. ( 5.16; ), satisfies the initial conditions of 
Eq. fl2.16D , namely, 

tf({-0, 7 )- 7 and !»E*&%1 (5.17) 



For large ( Eq. (5.15) leads to 



£ 2 



(5.18) 



which does not depend on u. Using this limiting behaviour we can suggest the following extrapolation 
formula 



£ = v 7 ^ / dC' 



C 



1 



1) 



Y J + 2 ln(C/ 7 ) 



(5.19) 



+ (exp(-C') 

= 2 + 2 ln(C/7) (5.20) 

In the first integral in Eq. (|5.19| ) the small values of does not contribute and we can replace it by the 
asymptotic al behaviour at large 

Using Eq. ( 5.16| ),Eq. ( |5.18| ) and Eq. ( |5.20D we can suggest an approximation for the solution of 
Eq. (|5.15|) , namely. 

C(£;7) = 7(l + e/2) + e/8 (5.21) 

Substituting Eq. (5.21) into Eq. ( |5 . 1 2| ) and choosing £' = £/2 that correspond to the smallest corrections 
(see ) , we obtain the final result 

N MPSI {0 = 1 _ exp (_ e 2 /16 _ 7(1 + e/2) ) (5 22) 

Therefore, the MPSI approximation leads to the geometrical scaling which satisfies the initial conditions 
at £ = 0. The £ 2 term does not depend of the scattering amplitude 7 oc a| while the next terms show such 
dependence. Comparing with the solution of Eq. (5.10) one can see that both have the same £ 2 dependence 
but the coefficients in front of of this term is quite different. We need to recall that the MPSI approach 
cannot be correct at very high energies |B3]. 
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6. Conclusions 



The main result of this paper is the reduction of the complicated problem of summing all Pomern loops to 
a Pomeron calculus in the zero transverse dimension but for 'time' y given by y = £ 2 /8 at £ S> 1. The 
principle features of this solution is the fact that the scattering amplitude depends on the only one variable: 
Q 2 (x)r 2 where Q s (x) is the saturation scale while r is the size of the scattering dipole. Such geometrical 
scaling behaviour has been proven for the mean field approximation and has been seen experimentally 
[ 31 1 . However, this scaling behaviour was not expected for the theory given by the functional integral 
of Eq. fll.ip . Indeed, at first sight, such theory can be rewritten as the Langevin equation for directed 
percolation ||3q| , namely, 

d _l^pXl = g| d *k>K(k,k')$(b,k';Y) - 2 -^.^(b,k;Y)^b,k;Y) + C(b,k;Y) (6.1) 
with 

< \((b,k;Y)\ > = 0; (6.2) 
< \((b,k;Y)({b',k';Y')\ > = l^.$(6,fc ; Y) (1 - $(b, k; Y)) 5 (2) {b-b') k 2 6® {k - k') 8 (Y - Y') 

1 » c 



As it has been shown in Refs. [22, ^J, 24] we rather expect from the statistical physics analogy a 
solution of the following type 

(~ \ 1 f , 9 9x , , 9 , 9sN ( (HQ 2 S R 2 )~ <HQ 2 R 2 ) >) 2 \ , N 

N^;Y) = ^ J dln(Q 2 R 2 ) N{ln(Q 2 /k 2 )) exp f- 1 Ws 2a 2 {Y) ) (6 ' 3) 

where N (ln(Q 2 /k 2 )^ is the scaling solution in the mean field approximation; < InQ 2 > is the average 
saturation scale which differs from Q s given by Eq. (|2.6| ) or Eq. (2.7), namely, f22]] 



< ln(Q 2 R 2 ) > = HQ 2 R 2 ) - ^ ) 7' Y (6.4) 



cr J 

~ ln 2 (V4) 



and variance of the saturation momentum a is proportional to 

a 2 =< ln 2 (Q 2 R 2 ) > - (< ln(Q 2 R 2 ) >f oc " s — Y (6.5) 

The Pomeron loops change drastically the geometrical scaling behaviour of the MFA solution in an 
obvious contradiction with our solution. Indeed, Eq. ( |6.3D means that the amplitude is a function of one 
variable 

Njr) ■■>'"'»- <M<S(r)«-)>'. (M) 



Despites that Eq. ( |6. 6|) looks very plausible we should not forget that it is correct for the functional 
integral of Eq. ( |1.2| ) which does not treat the impact parameter behaviour in consistent way. As has been 
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demonstrated in section 3 Eq. ([^), which is needed for a proof of Eq. (fb^), cannot be valid in the Pomeron 
loops. 

The reason for such difference between our approach and Eq. (|6.6| ) is clear from the consideration in 
section 3: we found that not all Pomeron loops lead to a violation of the geometrical scaling behaviour. 
Some of them contribute to a change of the scaling solution in the MFA. 

In our approach we have the following equation for the directed percolation instead of Eq. ( |6.l[) 

4^|^ = **(o - ^M.e*(o*(0 + ceo (6.7) 



with 



< |»n*C(OI > = 0; (6. 

47T<5? 



< ic(o c(Oi > = z j^f ■ z *co (i - - a 

These equations are equivalent to one dimensional Pomeron theory. The first and, perhaps, the only 
that this theory does, it leads to replacement of the MFA solution in Eq. ( |6.3[) by the solution that has 
been discussed here. Eq. fl6.8|) suggests that the correlation term in Eq. (|6.2| ) has a more complicated 



structure. We propose to sum Eq. ( |6.8| ) and Eq. (6.2) but , at the moment, we do not understand the 
correct b dependence in Eq. ( |6.2[ ). 

We would like to stress that we suggest a solution only for £ 3> 1 and matching with the initial 
conditions was done in an approximate approach (the MPSI approach). Therefore, in principle, could be 
that the solution we advocate here, just a wrong one but we believe that the MPSI approximation works 
well enough to be an argument against such a possibility. 

On the other hand, Eq. ( |6.3| ) is proven only in the vicinity of the saturation scale where, indeed, 
N (ln(Q2A 2 )) is tne solution for MFA which is averaged with a Gaussian weight. However, it is difficult to 
understand how Eq. ( |6.6| ) can be matched with the perturbative QCD solution of Eq. (|2.15| ) or Eq. fl2.16|) 
outside the saturation kinematic region. Indeed, din N/dln(k 2 R 2 ) = (l/a(Y)) dlnN/dr <C 1 instead of 
dlnN/dln(k 2 R 2 ) w 1. 



At first sight the numerical calculation in Refs. 15, 22, 23, 24, 37] can give us an information on how 



the solution behaves deeply inside of the saturation region. However, it is not the case at the moment 
, since: (i) all calculations have been performed for the diffusion kernel of Eq. (|2.5| ) (in other words, for 
the F-KPP equation) that cannot pretend to describe approaching to the asymptotic behaviour; (ii) the 
hope that such calculation can be generalized to the exact BFKL kernel is based on some conditions (see 
Rev. [37]]) which are plausible only in vicinity of the saturation scale; (iii) it has not been checked that 
the numerical solution reproduces the solution of Eq. ( p.!5| ) or Eq. ( p,16| ) , rather from fig. 16 of Rev. [37] 



one can conclude that they do not approach them; and (iv) the accuracy of the numerical calculations 
has not been checked by a comparison with the exactly solved model of the BFKL Pomeron calculus in 
zero transverse dimension. For the BFKL Pomeron calculus in zero transverse dimension the numerical 



calculation should coinside with the exact solution given in Ref. [29] with Eq. ( |6.2[) and should lead to 
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decreasing amplitude at high energy (see Ref.[38| and references therein) for Eq. ( |6^ ) without <3? 2 (6, k\Y) 
term. 

It should be stressed that the experimental data do not show any deviation from the geomerical scaling 
behaviour while the estimates with as ~ 0.25 show that we expect a considerable violation of this scaling 
behaviour accordingly to Eq. ( |6 . 3| ) . 

Summarising we do not think that there is any proof that our solution is inconsistent but , of course, 
we need to remember that this solution cannot be justified in the vicinity of the saturation scale. The only 
that we showed is that our solution does not contradict the exact solution in the perturbative QCD region 
(see Eq. (|1|) or Eq. (|2~16|) ). Such strate gy was successful in the case of the mean field approximation 
|p7[. On the other hand, solution of Eq. (|6.3|) is able to describe the system in the vicinity of the saturation 



scale but it is not hundred percents clear how this solution could be matched with Eq. (2.15) or Eq. ( |2.16|) . 
We would like to stress again that the most probable way out is that we need to replace the MFA solution 
by ours in Eq. (|fO|). 

We are working to include the exact BFKL kernel in our consideration and to understand better the 
kinematic region of matching with the perturbative QCD solution. This is a very important issue since the 



model kernel of Eq. (2^) does not give a correct high energy solution in the MFA (see Ref. f^). If we will 
use the correct solution in the MPSI approximation we will obtain the solution which being substituted in 
Eq. ( |6.3D leads to a geometrical scaling behaviour after integration over Xn(Q 2 s B?) . The only change is to 
replace Q 2 . in the definition of variable £ (see Eq. ( 2.1 1| ) ) by the new saturation momentum of Eq. ( |6.4[) . 

We hope that the paper will be useful in searching a solution in the saturation region beyond the 
scope of statistical physics analogy for the BFKL Pomeron calculus including the Pomeron loops. 
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